Choice of Integrator in the Hybrid Monte Carlo Algorithm
نویسنده
چکیده
We study efficiency of higher order integrator schemes for the hybrid Monte Carlo (HMC) algorithm. Numerical tests are performed for Quantum Chromo Dynamics (QCD) with two flavors of Wilson fermions. We compare 2nd, 4th and 6th order integrators at various quark masses. The performance depends on both volume and quark mass. On currently accessible large lattices ( V ∼ 24 ), higher order integrators can be more efficient than the 2nd order one only in heavy quark region, mqa > 0.3. Thus we conclude that for most full QCD simulations, except for heavy quark case, the usual 2nd order integrator is the best choice.
منابع مشابه
Higher Order Hybrid Monte Carlo at Finite Temperature
The standard hybrid Monte Carlo algorithm uses the second order integrator at the molecular dynamics step. This choice of the integrator is not always the best. We study the performance of the hybrid Monte Carlo algorithm for lattice QCD with higher order integrators in both zero and finite temperature phases and find that in the finite temperature phase the performance of the algorithm can be ...
متن کاملOptimal Scheduling of Battery Energy Storage System in Distribution Network Considering Uncertainties using hybrid Monte Carlo- Genetic Approach
This paper proposes a novel hybrid Monte Carlo simulation-genetic approach (MCS-GA) for optimal operation of a distribution network considering renewable energy generation systems (REGSs) and battery energy storage systems (BESSs). The aim of this paper is to design an optimal charging /discharging scheduling of BESSs so that the total daily profit of distribution company (Disco) can be maximiz...
متن کاملAlgorithms for Dynamical Fermions
This is the write-up of three lectures on algorithms for dynamical fermions that were given at the ILFTN workshop ‘Perspectives in Lattice QCD’ in Nara during November 2005. The first lecture is on the fundamentals of Markov Chain Monte Carlo methods and introduces the Hybrid Monte Carlo (HMC) algorithm and symplectic integrators; the second lecture covers topics in approximation theory and the...
متن کاملTesting and tuning new symplectic integrators for Hybrid Monte Carlo algorithm in lattice QCD
We examine a new 2nd order integrator recently found by Omelyan et al. The integration error of the new integrator measured in the root mean square of the energy difference, 〈∆H2〉1/2, is about 10 times smaller than that of the standard 2nd order leapfrog (2LF) integrator. As a result, the step size of the new integrator can be made about three times larger. Taking into account a factor 2 increa...
متن کاملAdaptive Step Size for Hybrid Monte Carlo Algorithm
We implement an adaptive step size method for the Hybrid Monte Carlo algorithm. The adaptive step size is given by solving a symmetric error equation. An integrator with such an adaptive step size is reversible. Although we observe appreciable variations of the step size, the overhead of the method exceeds its benefits. We propose an explanation for this phenomenon.
متن کامل